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ABSTRACT 

We present a comprehensive study to determine if the LINER/H II region transition spectrum 
in NGC 4569 can be generated solely by photoionization by the nuclear starburst. A review of 
the multiwavelength data from the literature reveals no additional sources that contribute to the 
ionization. We find that the young starburst dominating the UV emission is distinct from the 
nuclear population of A supergiants identified in the optical spectrum by Keel (1996). Spectral 
synthesis analysis provides constraints on the physical nature of the starburst, revealing a 5-6 Myr, 
approximately instantaneous starburst with subsolar metallicity. These results are used to model 
the spectral energy distribution of the ionizing continuum. Luminosity constraints place limits 
on the steepness of the extinction curve for the young starburst. The Savage & Mathis (1979) 
curve satisfies all luminosity constraints and the derived reddening is similar to the emission line 
reddening. These results imply extreme conditions in the nuclear starburst, with ^5xl0 4 O and 
B stars compacted in the inner 9" x 13"region of the nucleus. Using photoionization analysis and 
employing all observational constraints on the emission line gas, we find very specific conditions 
are required if the spectrum is generated solely by stellar photoionization. At least two spatially 
distinct components are required - a compact region with strong O III emission and an extended, 
low density component emitting most of the S II flux. A high density component is also needed to 
generate the O I flux. Additionally, a limited contribution from Wolf-Rayet stars to the ionizing 
SED is necessary, consistent with the results of Barth & Shields (2000). We present a physical 
interpretation for the multi-component emission line gas. 

Subject headings: galaxies: individual (NGC 4569) — galaxies: nuclei — galaxies: starburst 



1. Introduction 

One of the most basic properties of low ionization emission line galactic nuclei, their ionization and 
excitation mechanism, is not currently well understood. These objects include low ionization nuclear emission 
line regions (LINERs), which have optical spectra dominated by forbidden lines from low ionization species 
([O I] A6300, [S II] A6726, [N II] A6584), and "transition objects", with optical line ratios intermediate 
to H II regions (dominated by Balmcr lines and sometimes [O III] A5007) and LINERs. Various ionization 
mechanisms have been proposed for these nuclei, including photoionization by a low luminosity active galactic 
nucleus (AGN) (Ferland & Netzer 1983), stellar photoionization (Terlevich & Melnick 1985; Filippenko & 
Terlevich 1992; Barth & Shields 2000), and shock-heating (Dopita & Sutherland 1995), however, the true 
nature of only a few individual nuclei has been demonstrated conclusively (Barth et al. 1996; Ho et al. 1996; 
Gabel et al. 2000). 
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The physical nature of transition objects remains particularly uncertain. While many LINERs have 
broad emission line components (Ho et al. 1997b) and compact, flat radio cores (Nagar et al. 2000) indicating 
they are bona fide AGNs, transition objects do not generally exhibit these features. Thus, there is no direct 
evidence linking them to AGNs. If stellar photoionization is the dominant mechanism in transition objects, 
an explanation for the uncharacteristically strong low ionization lines is necessary. Recent studies have 
determined the essential features required for stars to generate such a spectrum: (a) The spectral energy 
distribution of the stars must be "hard" compared to typical H II regions, i.e., it must have stronger relative 
emission beyond the He I and He II ionization limits. A harder SED produces an extended partially ionized 
zone in the emission line region, thereby enhancing the LINER features (Pequignot 1984; Shields 1992; 
Barth & Shields 2000). (b) Also, a low ionization parameter is required (Shields 1992; Barth & Shields 
2000), which is a function of the starburst's luminosity, nebular density, and the geometry of the nucleus. 

The transition object nucleus in NGC 4569 is an ideal object for testing a possible link between starburst 
activity and low ionization nebular emission. Of all LINERs and transition nuclei with available UV data, 
NGC 4569 exhibits by far the most pronounced nuclear starburst activity. Maoz et al. (1998) found that 
the ionizing luminosity of the nuclear starburst in NGC 4569 is sufficient to generate the observed emission 
line luminosity. Barth & Shields (2000) demonstrated that a large population of very hot Wolf-Rayet stars 
could harden the integrated SED of a starburst enough to produce the partially ionized zone required to 
produce a transition object spectrum. They also found that the estimated age of the stellar population in 
NGC 4569 roughly coincides with the W-R phase of evolution in a starburst. 

However, there are no observational signatures of W-R stars in the spectrum of NGC 4569 (Barth & 
Shields 2000). Additionally, W-R galaxies typically have emission line spectra characteristic of H II regions, 
rather than transition objects or LINERs (Guseva et al. 2000). Furthermore, it is unclear if the ionization 
parameter in NGC 4569 is sufficiently low to generate the observed spectrum. Hence, the questions remain: 
Is the observed line emission in NGC 4569 generated solely by photoionization by the nuclear starburst? If 
so, what physical conditions in the nebular gas and starburst are required? 

We have undertaken a detailed analysis of the emission in NGC 4569 to answer these questions, incorpo- 
rating all available observational constraints. In § 2, we summarize the multiwavelength data of NGC 4569 
from the literature and present the archival data used in our analysis. In § 3, we apply population syn- 
thesis analysis to these data to place constraints on the nature of the nuclear starburst and the emitted 
flux distribution. Using the derived constraints on the SED and luminosity of the starburst, along with 
observational constraints on the emission line gas, we use photoionization modeling in § 4 to determine if 
stellar photoionization can generate the observed transition object spectrum emitted by NGC 4569. Finally, 
we discuss the implications of our results in § 5. 

2. Nuclear Data and Multiwavelength SED in NGC 4569 

In the following, we describe the optical and UV data used in our analysis and summarize results 
obtained from the literature for several bandpasses. These multiwavelength observations are crucial for 
determining the physical mechanisms that power the nuclear emission. 

For our study, we adopt a distance of 16.8 Mpc to NGC 4569, given by de Vaucouleurs et al. (1991). Due 
to the anomalous radial velocity of NGC 4569, v — —223 km s _1 , there has been much debate concerning 
this distance, and it remains highly uncertain (Stauffer et al. 1986, and references therein). We explore the 
affects of this uncertainty on our results in § 5. 
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2.1. Spectrum and Morphology of the Emission Line Region: A LINER/H II Region 

Transition Object 

NGC 4569 was observed as part of the Palomar survey of nearby galactic nuclei by Ho et al. (1997a) 
using the Hale 5m telescope with a 2" x 4" aperture («160 x 320 pc). Based on the observed emission line 
ratios, [O I] A6300/Ha, [S II] A6726/Ha, [N II] A6584/Ha, and [O III] A5007/H/?, they classified NGC 4569 as 
a LINER/H II region transition object nucleus. The observed line luminosities and estimated uncertainties 
measured by Ho et al. are given in Table 1. Using the Savage & Mathis (1979) extinction curve and assuming 
an Ha/H/3 flux ratio of 2.9, which is typical of emission line nuclei (Osterbrock 1989), we derive a reddening 
of E(B — V) =0.50 (±0.15) for the emission lines. Dereddened line luminosities are presented in column 3 
of Table 1. The quoted uncertainties were derived from the uncertainties in the flux measurements and the 
reddening, using propagation of errors. The uncertainty in the extinction curve in NGC 4569 introduces an 
additional uncertainty and is examined in more detail in § 5.2. 

The morphology of the nuclear line emission in NGC 4569 is delineated in the Hubble Space Telescope 
(HST) WFPC2 narrowband imagery presented by Pogge et al. (2000). The continuum subtracted Ha 
(F656N) and [O III] A5007 (F502N) images shown in Figures lc and 3b of their study reveal a bright 
nuclear core with total extent < 1" (80 pc). There is also a more extended, diffuse nebular emission 
component apparent in both images. Examination of the spectral bandpass of the F656N filter reveals little 
contamination from the [N II] AA6548,6584 lines at the blucshift of NGC 4569, and is therefore a good measure 
of the Ha emission. For each image, Pogge et al. have measured the flux emitted in the inner R =0".23 (20 
pc) core of the nucleus, as well as the more extended nuclear region, R =2". 5 (250 pc). Their results give 
a measure of the compactness of the line emission: L2o pc / ^250pc =0.31 for Ha and L20pc/L25o pc =0.74 for 
[O III]. These flux ratios provide valuable constraints for our photoionization analysis in § 4. 

2.2. UV Data: Evidence for a Compact Nuclear Starburst in NGC 4569 

We obtained Faint Object Spectrograph (FOS) / HSTW data from the Space Telescope Science Institute 
(STScI) data archives. These data were originally published by Maoz et al. (1998). Two exposures with the 
G130H grating, and one each with the G190H and G270H gratings were used. These data have a spectral 
resolution of AA ^1-2 A . Details of the observations are summarized in Table 2. The two G130H grating 
exposures were coadded by weighting them according to their exposure times. A combined spectrum was 
then created from the three spectra and corrected for the systemic velocity of the galaxy. 

All exposures were taken through the 0".86 square aperture of the FOS. In a study of ffi>T/WFPC2 
images taken with the F218W filter (A centrai =2200 A), Barth et al. (1998) found the UV nucleus of NGC 4569 
to be very compact though clearly resolved, with a FWHM of 0".16 x 0".ll (13 x 9 pc 2 ). Thus, the FOS 
data sample the entire compact UV emitting region. Furthermore, we find that the UV flux of the compact 
nucleus measured from the HST imagery by Barth et al. matches the International Ultraviolet Explorer 
(IUE) spectrum obtained through the 10" x 20" aperture (Keel 1996) within the noise of the data. Hence, 
our comparison shows no significant UV emitting sources, other than the compact nuclear starburst, are 
contributing to the ionization of the emission line gas. 

Maoz et al. (1998) found that the FOS UV spectrum of NGC 4569 is dominated by a very young 
starburst. The resonance doublets of C IV AA1548,1551, N V AA1239,1243, and Si IV AA1393,1403 exhibit 
P Cygni profiles, indicative of the strong stellar winds produced in O and B stars. Maoz et al. found that 
the line profiles of the wind features in NGC 4569 are similar to the starburst in NGC 1741, which has a 



-4 - 



derived age of ~5 Myr (Conti et al. 1996). They also identified several photospheric features that appear in 
young, hot stars. 

Numerous interstellar lines, both Galactic and intrinsic to the host galaxy, are also apparent in the 
UV spectrum of NGC 4569 (Maoz et al. 1998). Strong features from the multiplets of C I near A1277, 
A1280, A1329, and A1561 were detected at the systemic velocity of NGC 4569, which implies the starburst 
spectrum is seen through a large column of cool, possibly high density gas intrinsic to NGC 4569. Thus, 
large extinction of the UV spectrum due to dust is expected. Further evidence for a high extinction of the 
starburst comes from the flatness of the observed UV spectrum between 1200-3000 A (Maoz et al.), since an 
unreddened population of O and B stars would have a sharply rising continuum toward shorter wavelengths. 

Besides the red-shifted emission lobes associated with the stellar P Cygni resonance profiles, there are 
no other apparent emission lines in the UV FOS spectrum of NGC 4569 in the 1200-3000 A range. Notably, 
there are no detectable stellar emission features that might indicate the presence of Wolf-Rayet stars, such 
as He II A1640 or C III] A1909, nor are there UV nebular emission lines. Upper limits on nebular emission 
line fluxes were determined by Maoz et al. (1998) and are included with our optical fluxes in Table 1. 



2.3. The Optical Counterpart to the UV Starburst 

The continuum flux of the optical emission coinciding with the UV starburst can be derived from analysis 
of HST/WFPC2 continuum images. Combined with the UV spectrum, this provides constraints in deriving 
the extinction of the central starburst (see Section 3.3). Exposures taken with the F547M (A centra / =5470 A) 
and F791W (A centra/ =7920 A) filters were used because their bandpasses are sufficiently narrow to exclude 
contribution from strong nebular line emission, thereby providing fairly reliable stellar continuum fluxes. 
Although the total CCD detector counts are high, data for the bright core in each of these images remain 
unsaturated. We integrated the total fluxes in these images over a 0".35 x 0".35 spatial region corresponding 
to the approximate size of the UV nucleus as determined from the WFPC2 F218W image (Barth et al. 1998). 
In Figure 1, we show the F547M image and delineate the region on the sky used for our core flux measurement. 
Specific fluxes were derived using the conversions provided by Holtzman et al. (1995), assuming a flat flux 
distribution. This yields F ° =5.0xl0~ 15 ergs s -1 cm~ 2 A -1 and F ° =2.6xl0~ 15 ergs s _1 cm~ 2 A -1 . 

J 5470A ° 7920A b 

In a ground-based study of NGC 4569 using a 2" x 4" aperture with the 2.5m Isaac Newton Telescope, 
Keel (1996) deduced that the optical stellar continuum is dominated by a population of A-type supcrgiants 
from the strong narrow Balmer absorption features. These stars point to a much older population (> 15 Myr) 
than implied by the UV data. Figure 1 reveals that the optical stellar continuum emission is clearly more 
extended than the compact UV starburst. We have measured the flux in the inner 2" x 4" region of the 
F547M image and find that <20% of the optical flux coincides with the UV source, i.e., the inner 0".35 x 0".35 
region measured above. Therefore, the population of A supcrgiants must be in the extended component seen 
in Figure 1 and cannot be directly associated with the compact starburst. This indicates that the nucleus 
of NGC 4569 contains at least two distinct young stellar populations: (1) the very young, very compact UV 
core and (2) a spatially more extended population dominated by A supergiants, with few, if any, young O 
and B stars present. 
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2.4. Multiwavelength Observations 

Although the UV spectral data and imagery imply that a very young stellar population is dominating 
the emission in NGC 4569, this does not rule out the possibility that an obscured AGN is present, which 
can contribute to the ionization in the emission line gas. Below, we summarize multiwavelength observations 
obtained from the literature to test this possibility. 

Radio Emission 

There are no signatures of an AGN in radio observations of NGC 4569. High angular resolution radio 
data obtained with the Very Large Array (VLA) at 6 cm and 20 cm reveal only weak, extended («4", 
320 pc) emission in the nuclear region, with integrated fluxes of F§ cm =1.5 mJy and Fi^ cm =10 mJy (Neff 
& Hutchings 1992; Hummel et al. 1987). No compact core component is detected in this data, with an 
upper limit on any unresolved flux of FQ Cm <0.4 mJy at 1" angular resolution (Neff & Hutchings 1992). 
The spectral index of 0:20-6 = —1-5, derived by extrapolating between the fluxes at A =20 and 6 cm and 
assuming F v oc v~ a , is much steeper than the slope expected for pure free- free emission in an H II region, 
ctff, revealing a substantial non-thermal contribution. We estimated the contribution of free- free emission to 
the observed flux using the Ha emission line luminosity given in Table 1 , following the method of Osterbrock 
(1989). Assuming all of the Ha emission is generated by photoionization in a nebula at T e =10 4 K gives 
^0.7 mJy for the free-free flux, or roughly one-half of the observed value at 6 cm. Subtracting this from the 
observed fluxes implies a spectral index of 020-6 = —1-9 for the non-thermal emitting component. 

The steep slope of the radio emission in NGC 4569 may be the signature of supernovae remnants (SNRs). 
For example, a survey of SNRs in M31 shows a range of spectral slopes, a 2 o-6 ~ —0.35 to —1.6 (Dickel & 
D'Odorico 1984). Assuming the non-thermal component of radio emission in NGC 4569 is from SNRs and 
using an SNR luminosity at 6 cm of L 6cm =l.lxl0 24 ergs s _1 Hz -1 (Dickel & D'Odorico ), approximately 
100-200 SNRs are required in the nucleus of NGC 4569, depending on the contribution of the free- free thermal 
emission. 

Far-Infrared Emission 

Far-infrared (FIR) fluxes at A = 12, 25, 60, and 100 /xm are available from IRAS observations (Neff & 
Hutchings 1992). These data, at ~30" resolution, sample an extended region in NGC 4569. Since there 
are no signatures of a central AGN in the radio data, the FIR emission most likely arises from warm dust. 
Assuming most of the nuclear emission is reprocessed as thermal dust emission, the total FIR luminosity, 
Lfir =2x10 ergs s^ 1 , provides an upper limit on the bolometric luminosity of the starburst. We have 
compared the FIR spectral energy distribution (SED) in NGC 4569 to blackbody flux distributions at various 
temperatures and find that at least two dust temperatures are required to fit the data. This is not surprising 
considering the large spatial region sampled by the IRAS data. The warmer dust may be associated with the 
inner nuclear region, close to the energetic starburst and possibly in the same gas giving rise to the observed 
C I absorption, while the cooler dust may be located in more extended regions. 

X-Ray Emission 

Previous studies of the nuclear X-ray emission in NGC 4569 also reveal no signatures of an AGN. 
Terashima et al. (2000) found no compact hard X-ray emission in ASCA observations, placing an upper limit 
of L2-10 keV <3xl0 39 ergs s _1 on any unresolved core luminosity. Thus, if an AGN capable of powering the 
nebular line emission in NGC 4569 is present, it would have to be heavily obscured. A compact component 
has been detected at soft X-ray energies with ROSAT/HRI, Lo.1-2.4 keV =6xl0 39 ergs s^ 1 (Colbert & 



- 6 - 



Mushotzky 1999), and with ASCA, L . 5-2.0 keV =6xl0 39 ergs s _1 (Terashima et al.). The ~4" resolution of 
ROSAT/RRl imposes the most rigid constraint on the size of any compact emitting region, but still samples 
the entire spatial region of the compact stellar core, the nebular gas, and the radio emitting source. Spectral 
analysis of the soft ASCA data by Terashima et al. indicates that it is not consistent with scattering of non- 
thermal emission by an obscured AGN. Combining this result with the low upper limit on the compact hard 
X-ray flux, they concluded that the unresolved soft X-ray emission is associated with the nuclear starburst 
rather than an AGN. 

Young starbursts are known to emit soft X-rays. The combined emission from SNRs, interstellar gas 
heated by supernovae and stellar winds, and X-ray binaries can all contribute to the integrated X-ray 
emission in a starburst (Fabbiano 1989, and references therein). Unfortunately, distinguishing between 
these mechanisms in NGC 4569 by spectral modeling is not possible with the available data. The core 
luminosity in NGC 4569 is consistent with other starburst nuclei, which are typically in the 10 39 ergs s _1 
range (Fabbiano ). However, the stellar population and the X-ray source are more compact in NGC 4569 
than in typical starburst nuclei. 

The Broadband Spectral Energy Distribution 

In Figure 2, the observed SED of the nuclear emission in NGC 4569 from radio to soft X-ray energies is 
plotted. UV fluxes were obtained from the FOS spectrum. The optical data are for the inner 0".35 x 0".35 
core of emission derived in § 2.3, while the radio data points are upper limits to any unresolved flux. Due 
to the low angular resolution of IRAS, the FIR data represent upper limits to the nuclear core flux. The 
soft X-ray SED in Figure 2 was derived by fitting a power law flux distribution (F v oc v~ a = —1.5) to the 
observed compact (within «4") ROSAT/BRl flux quoted in Colbert & Mushotzky (1999). This was done 
only as an approximation and has no physical basis. 

Qualitatively, the entire broadband nuclear SED is consistent with emission associated with a young 
stellar population. Figure 2 shows that the emitted power is strongest in the optical-UV region, and possibly 
in the FIR (depending on the contribution from the inner nucleus), with relatively low luminosities at radio 
and X-ray frequencies. The spectral signatures in the UV indicate a nuclear starburst is dominating the 
emission and the FIR luminosity can be understood solely as reprocessed stellar radiation with no indications 
that an AGN is present. 

3. Spectral Population Synthesis of the Nuclear Starburst 

To test if photoionization by the nuclear starburst can generate the LINER/H II region transition 
spectrum emitted in NGC 4569, we must first establish the integrated spectral energy distribution and 
luminosity of the stars. Thus, we employ spectral population synthesis analysis in this section to derive 
constraints on the evolutionary stage and physical conditions in the starburst. These results are then 
combined with stellar atmosphere models to generate a model SED for subsequent photoionization analysis. 
The derived flux distribution, along with observational constraints described in the previous section, place 
limits on the reddening, and hence the luminosity of the nuclear starburst. 
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3.1. Description of the Modeling Code 

In the following, we employ a population synthesis modeling code, CLUST_FLUX, which generates 
UV spectra for a starburst of specified age, metallicity, and initial mass function (IMF). CLUST_FLUX 
computes the present day stellar population using the stellar evolutionary models of Schaller ct al. (1992), 
Schaerer et al. (1993a, b), Charbonnel et al. (1993), and Meynet et al. (1994). These models predict the 
isochrones in luminosity - effective temperature parameter space (L — T e //) for a given age and metallicity. 
Combining these results with the IMF, the relative number of stars for all locations in the L — T e ff plane is 
then fully specified. Stellar evolutionary tracks that include the effects of mass loss are available for various 
metallicitics, ranging from 0.05-2 Z Q , where Z Q represents solar metallicity. 

CLUST_FLUX generates a composite model starburst spectrum using short wavelength IUE data (1170- 
2000 A) . A weighted sum of individual stellar spectra is created from a library of IUE data compiled for the 
Milky Way and Large and Small Magellanic Clouds, which are classified according to effective temperature 
(T e ff) and surface gravity (log(g')). Thus, this library is a 3-dimensional grid of spectral data in T e ff, 
log(g), and Z, covering three metallicities: Z Q , 0.25 Z Q (LMC), and 0.1 Z (SMC). For stellar types not 
corresponding to exact matches in T e ff — log(g) — Z in the IUE library, interpolated spectra are generated. 
To correct for extinction, with often unknown, highly variable reddening laws for stars in the IUE data, all 
spectra have been matched to the flux distributions of Kurucz (1992) model atmospheres for corresponding 

T eff - log(g) - Z. 

The UV spectral classification criteria of Smith Neubig & Bruhweilcr (1997, 1999), aided with the UV 
line identifications of Dean & Bruhweilcr (1985) and Bruhweiler et al. (1981), were used to classify the 
UV stellar spectra for OB stars in the SMC and LMC. In practically all cases, these classifications were in 
close accord with corresponding visual spectral types in all three galaxies. High resolution IUE spectra of 
Galactic stars are available for the solar metallicity grid and have been binned into AA =0.25 A intervals. 
Since there are few high dispersion IUE spectra for stars in the Magellanic Clouds, the stellar library of 
sub-solar metallicity SMC and LMC stars is comprised of low resolution (AA~7 A ) IUE spectra. 

Due to the large uncertainties involved in predicting the population of Wolf-Rayet stars in a star- 
burst (Schmutz & Vacca 1999), we do not include them in the integrated model starburst spectra. Thus, 
CLUST.FLUX does not predict the strengths of the broad UV emission features that are signatures of 
Wolf-Rayet stars, such as He II A1640, C IV A1550, and C III] A1909. However, emission from W-R stars 
can significantly alter the integrated ionizing emission in a starburst and has important consequences for 
photoionization analysis, as demonstrated by Barth & Shields (2000). Therefore, we explore observational 
constraints on the W-R population in NGC 4569 separately in § 3.2. 

3.2. Stellar Population Synthesis Results 

Our method is to examine the effect of age on an instantaneous burst of star formation using standard 
model parameters and to determine the best fit to the spectrum of NGC 4569. The effects of varying 
metallicity, the IMF, and star formation rate can then be explored to better constrain the physical nature 
of the starburst. In the following analysis, we adopt the high mass-loss rates of the Meynet et al. (1994) 
evolutionary models. 

Age Constraints 

We have generated model spectra for an instantaneous burst of star formation for a range of ages using 
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the standard IMF of Miller & Scalo (1979) and solar mctallicity. Their IMF, £ms, is defined by a three 
component power law function of mass, with spectral indices a\ = —0.4, a 2 = —1.5, and a 3 = —2.3 in the 
mass ranges Mi =0.1-1 M©, M 2 =1-10 M Q , and M 3 =10-120 M Q , respectively. 

The synthesized, normalized spectra in four far-UV spectral regions for instantaneous burst models with 
ages 0, 2, 4, 5, 6, and 8 Myr are plotted in Figure 3. These show the spectra for the C IV, N V, Si IV, and 
He II features, as well as important photospheric diagnostic lines. The normalized spectrum of NGC 4569 
is superposed on all model spectra. Due to the limited S/N in the FOS spectrum of NGC 4569, we have 
convolved the data with a a =1.0 A Gaussian filter for clarity. To reproduce the spectral resolution of the 
smoothed FOS data, all model spectra were smoothed using a Gaussian filter with a =1.5 A, which reflects 
the best match between the spectral widths of the photospheric features in the model and observed spectra. 
The uncertainties in our diagnostic line fits are dominated by uncertainties introduced by the photon noise 
in the spectrum of NGC 4569. 

In Figures 3a-c, the stellar wind profiles of the N V, C IV, and Si IV resonance doublets are seen to 
vary dramatically with age between and 8 Myr. The observed N V AA1239,1243 feature in the spectrum 
of NGC 4569 does not exhibit the broad blue absorption wing seen in the very young starburst models, 
indicating an age of at least 4 Myr. However, the red-shifted N V emission lobe of the P Cygni profile is 
clearly present in NGC 4569 and is best reproduced by the 5-6 Myr instantaneous burst models. Since the 
Si IV AA1394,1403 doublet feature is most pronounced in high luminosity supergiants, it is a very useful 
age discriminate. Strong Si IV appears for only a short range of ages and does not occur in the youngest 
starbursts, as illustrated in Figure 3b. The best-fit to the Si IV wind feature is 4-5 Myr, with a rigid upper 
limit of 6 Myr. For later ages, the Si IV feature vanishes entirely. The C IV AA1548,1551 P Cygni feature 
is strongest in the very youngest starbursts and weakens systematically with age (Figure 3c). It is best 
matched by the 5 and 6 Myr models. The subordinate line, N IV A1718, which exhibits a P Cygni profile 
in spectral types 07 and earlier, is not present in NGC 4569. Our models predict this feature will become 
undetectable in starbursts older than 4 Myr. In Figure 3d, the spectrum for He II A1640 is plotted. Since the 
CLUST_FLUX models do not include a contribution from W-R stars, the model predictions for this spectral 
region are similar for different ages. No evidence for a broad He II A1640 emission feature is evident in the 
spectrum of NGC 4569, thereby limiting the number of W-R stars in the starburst (see discussion below). 

The strengths of the photospheric lines also provide important age diagnostics. The Si III AA1295-1303 
lines from the 3 P- 3 p° multiplct peak in strength in late B-type stars and, thus, increase through 8 Myr. 
The Si III lines that are uncontaminated by interstellar O I A1302 and Si II A1304 absorption are plotted 
in Figure 3a. In the NGC 4569 spectrum, the strength and width of the Si III A1294.5 multiplct line 
is inconsistent with the Si III A1296.7 and Si III A1298.9 lines and is therefore an unreliable diagnostic. 
However, the Si III A1296.7 and Si III A1298.9 line strengths are matched by models with ages > 5 Myr. The 
Si III A1417 line, shown in Figure 3b, exhibits the same trend with age as the Si III AA1295-1304 multiplct 
and indicates an age of 5 Myr or older in NGC 4569. In Figures 3a and b, the C III A1247 and C III AA1426- 
1428 line strengths also vary with age. The strength of each of these lines in NGC 4569 implies an age 
> 4 Myr. At later ages, the lines show little variation and provide no upper limit for the ages of interest 
here. In Figure 4, we plot the model and observed spectra for 1890-1990 AThis samples numerous Fe III lines 
(marked on the top of Figure 4) that become prominent in B supergiants, thus providing an additional age 
diagnostic. The Fe III multiplets are seen to increase in strength with later ages. The spectrum of NGC 4569 
is best matched at these wavelengths by ages of 5 Myr and older. 

The combined fits to the stellar wind and photospheric features indicate an age of 5-6 Myr for the 
starburst. Particularly strong limits come from the presence of the N V wind feature, indicating the starburst 
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is < 6 Myr, and the strength of the short-lived Si IV wind line, implying an age of 4-6 Myr. The photosphcric 
line diagnostics imply an age of at least 5 or 6 Myr for NGC 4569. 

Initial Mass Function 

Model spectra were also generated using different values for the IMF, star-formation rate, and metallicity 
in an attempt to constrain these parameters in NGC 4569. To test the IMF, we produced models using both 
a natter (a = —1.5) and steeper (a = —3) IMF relative to £ MS for the stellar masses M >10 M . Like 
the previous models, these spectra were generated for an instantaneous burst of star formation using solar 
metallicity. 

A flatter IMF results in a spectrum that is more heavily weighted with features from the most massive 
main-sequence (MS) stars. Thus, features that arise primarily in the most evolved stars are enhanced, 
resulting in stronger P Cygni features at all ages relative to the models using £ms- Conversely, we find the 
predicted strengths of the photospheric lines appearing in the spectrum of NGC 4569 are diminished for 
ages younger than 6 Myr in the flat IMF models. Thus, if the IMF in NGC 4569 is flatter than £ms, the 
observed strengths of the photospheric features indicate an earlier age than derived previously. 

Of course, the models having a steep IMF give opposite results. This provides a constraint on the IMF 
in NGC 4569. Specifically, since the observed strength of the Si IV wind feature in NGC 4569 approaches 
the maximum strength predicted by the standard models (Figure 3b), the IMF in NGC 4569 must be at 
least as flat as £ms- Otherwise, the Si IV feature would be underestimated in models for all ages. 

Star Formation History 

The star formation history of the starburst in NGC 4569 is also constrained by the observed strength 
of the Si IV wind feature. Since a pronounced Si IV wind feature occurs over only a short age interval and 
since the strength of this feature in NGC 4569 equals the maximum predicted strength for an instantaneous 
burst, the duration of star formation must be relatively short. For star formation lasting longer than about 
3 Myr, the Si IV wind feature would be diluted and weaker than observed in NGC 4569. A similar argument 
applies in limiting the contribution from multiple bursts of star formation. Thus, the starburst in NGC 4569 
is consistent with a single burst of approximately instantaneous star formation. Barth & Shields (2000) 
reached a similar conclusion based on the small spatial extent of the starburst. 

Metallicity 

Evolutionary models show that the mass-loss rates of massive stars are highly dependent on metallicity. 
Since mass-loss rates determine the stellar evolutionary tracks in the HR diagram and the rate of progres- 
sion along those tracks (Schallcr et al. 1992), the constituents of a stellar population are very sensitive to 
metallicity. 

With the LMC and SMC data in our stellar library, it is possible to test the effect of metallicity 
on a starburst spectrum. We have generated spectra for both Z =Z©/4 (LMC) and Z =Z Q /10 (SMC) 
for comparison with our standard models. The most striking difference is the decrease in strength of the 
wind features with decreasing metallicity. For example, the 5 Myr, Zq/10 model spectrum significantly 
underestimates the strengths of all three wind features observed in NGC 4569. Although the observed 
strengths of the N V and C IV features can be reproduced by Z =Z /1O models with earlier ages, the Si IV 
line cannot. Additionally, the weakness of the photospheric features in the low metallicity, young burst 
models are incompatible with the spectrum of NGC 4569. This implies Z >Z Q /10 for NGC 4569. Due to 
insufficient S/N in the spectrum of NGC 4569, we cannot place more rigid constraints on the metallicity. 
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Wolf-Rayet Stars 

The burst age of 5-6 Myr derived from our population synthesis analysis coincides with an age where 
large numbers of W-R stars may be present, although this depends strongly on the precise age and metallicity 
of the starburst (Meynet et al. 1994; Schmutz & Vacca 1999; Guseva et al. 2000). As discussed by Barth 
& Shields (2000), no W-R emission features are detectable in the optical spectrum obtained by Ho et 
al. (1997a). Furthermore, we find no broad He II A1640 (W\ <1.0 A) or C III] A1909 emission, the UV 
signatures of W-R stars, in the FOS spectrum of NGC 4569. Additionally, the C IV A1550 wind feature was 
well matched by our population synthesis models using only normal O and B stars (Figure 3c), which limits 
the contribution to this feature from W-R stars. 

Since the W-R population in a starburst is extremely sensitive to age and metallicity, the absence 
of detectable W-R features in NGC 4569 provides further constraints on the starburst. We employ the 
population synthesis results of Schacrcr & Vacca (1998), who derive the strength of the He II A1640 emission 
feature in an instantaneous starburst as a function of age and metallicity. Their models employ the same 
stellar evolutionary tracks used in CLUST_FLUX and incorporate a full treatment of the relative numbers of 
the different classes of W-R stars. In figure 10 in Schaerer & Vacca , they give the predicted equivalent width 
of the W-R He II A1640 emission line as a function of starburst age for various metallicities. This figure 
predicts W\~6 A for a 5-6 Myr instantaneous starburst with solar metallicity, which is six times greater 
than our measured upper limit. The He II emission is systematically reduced for lower metallicities, with 
Wa~1 A in the 5 Myr Z =Z Q /2.5 model, and smaller for older starbursts. Thus, if the spectral synthesis 
model predictions of the W-R population is correct, this indicates Z <Z /2.5 in NGC 4569. 

The SED, Luminosity, and Mass of the Nuclear Starburst 

Our spectral synthesis analysis indicates that a 5-6 Myr instantaneous burst of star formation with 
Z Q /10< Z <Z Q /2.5 and a Miller & Scalo (1979) IMF matches the overall spectrum of NGC 4569 quite 
well. CLUST_FLUX is similar to the well established population synthesis code of Leitherer et al. (1999), 
STARBURST99. We have computed numerous comparison models and find the spectra generated by the 
two models compare favorably. Specifically, they give identical results for the best fit to the starburst in 
NGC 4569. These results can now be used to construct a model of the ionizing flux distribution for the 
starburst using STARBURST99. Their models employ the model atmospheres of Lejeune et al. (1997), 
based on the plane-parallel, LTE, line-blanketing models of Kurucz (1992), for the OB stars and those of 
Schmutz et al. (1992) for the W-R stars. 

In Figure 5, we plot the 5 and 6 Myr model SEDs from STARBURST99 computed for Z =Z /2.5, con- 
sistent with our spectral synthesis results. A dramatic change in the shape of the SED occurs at wavelengths 
shortward of the He + ionization edge due to the disappearance of most of the hot W-R stars shortly after 

5 Myrs. 

As discussed in § 2, the strong interstellar C I and the flat observed UV continuum in the FOS spectrum 
implies a substantial extinction of the nuclear starburst emission. We adopted the extinction curve of Savage 

6 Mathis (1979) to correct for this. We address limits on the extinction law in NGC 4569 in more detail in 
§ 5.2. Comparison to the shape of the model SED between 1200-1600 A implies a reddening of E(£> — V r )=0.5 
of the starburst emission. This value is consistent with the reddening derived for the emission lines from 
the Balmer decrement. Integrating the model flux distribution then gives Lb i =1.5x10 43 ergs s _1 and 
L>ii?yd =10 42 ergs s _1 for the bolometric and ionizing luminosities, respectively. 

The nuclear starburst in NGC 4569 is quite remarkable. Integrating the IMF of Miller & Scalo (1979) 
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between 0.1 M Q <M* < 120 M© gives a total mass of 3xl0 7 M Q in the stellar population. There are ~5xl0 4 
luminous O & B stars (assumed to be those with MS masses M* >15 M©). Thus, for the massive stars 
alone, the cluster has a stellar density of about 10 pc~ 3 . We note that these results are sensitive to the 
uncertain extinction law of the starburst (see discussion in § 5.2). Comparison to the starburst nucleus 
in NGC 7714 (Gonzalez Delgado ct al. 1999) reveals NGC 4569 has about the same number of OB stars 
compacted in a volume that is 3000 times smaller. Though extreme, the stellar densities in NGC 4569 are 
not unprecedented. For example, the stellar density is believed to be similar in the Galactic center (Genzel & 
Townes 1987, and references therein). Furthermore, the densities and luminosities in the super star clusters 
in the Irregular galaxies NGC 1569 and NGC 1705 (O'Connell et al. 1994) are comparable. Interestingly, 
these latter objects are considered to be young globular clusters. This may indicate that the starburst in 
NGC 4569 is a gravitationally bound system. 

4. Stellar Photoionization Modeling of the Transition Nucleus in NGC 4569 

The primary question we are addressing is: Can ionization by a pure starburst explain the LINER/H II 
region transition object spectrum in NGC 4569? Now that we have a model of the SED and luminosity of 
the ionizing emission of the starburst, we can explore this using photoionization analysis. In the following, 
we present models using our derived SED and all the available observational constraints on the physical 
conditions of the emission line gas. 

4.1. Modeling Method, Assumptions, and Parameters 

Models were generated using the photoionization code CLOUDY (Ferland 1996). In computing these 
models, we assumed the stellar emission comes from a centrally located point source. Considering the simple 
geometry of the compact starburst evident in the WFPC2 UV image, namely a single, compact knot (Barth 
et al. 1998), this represents a good approximation to the actual morphology. 

Our method in this analysis was to find the simplest model that best reproduces the emission line 
spectrum observed in NGC 4569. Thus, we sought the fewest number of emission line gas components 
that combine to match the spectrum. Each model gas "component" was specified by a constant hydrogen 
gas density, n#, and distance from the ionizing starburst, R. Combined with the stellar luminosity, these 
parameters determine the ionization parameter for each component, 

U = 4:irL/n H R 2 c, 

which determines the ionization state of the gas. 

Modeling the emission line region as one or a few discrete components is a simplification. The distribu- 
tion of emission line clouds is certainly more complex than forming a set of distinct, perfect shells centered on 
the starburst. However, the minimal component method is the most instructive, since it reveals the general 
physical conditions that are required to generate the observed emission line spectrum. 

All models were computed for an open geometry, i.e., all emission escapes the clouds without interacting 
with emission line gas on the opposite side of the starburst. This assumes a priori a covering factor less than 
unity. Furthermore, all emission line clouds were assumed to be matter bounded, with an arbitrarily large 
column of neutral gas in the backside of the clouds. 
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The WFPC2 [O III] and Ha flux measurements by Pogge et al. (2000) provide constraints on the 
distribution of the emission line gas. They found that ~30% of Ra and ^75% of the [O III] flux is emitted 
in the inner R=20 pc region of the nucleus (see § 2.1). The electron density of the [S II] emitting region can 
be determined from the ratio of the [S II] A6716 to [S II] A6731 line fluxes (Osterbrock 1989). In NGC 4569, 
F[S//1 \67w/^[Sii] A673i = l-0 = ' = 0.2, which indicates an electron density of n e =600 +3OO _2oo cm" 3 . This is well 
below the critical density of these lines, n c ^10 4 cm" 3 , and therefore provides an important diagnostic for 
photoionization analysis. 

4.2. Model Results 

We computed a series of models with the 5 Myr starburst SED (Figure 5) at R=20 pc, consistent with 
the size of the [O III] A5007 emission core derived from the WFPC2 imagery. Guided by our population 
synthesis results, we assumed Z=Z Q /2.5 for the nebular gas, where Z Q is from (Grevesse & Anders 1989), 
modified by depletions onto dust grains. Gas-phase depletions and dust grain abundances typical of the 
Galactic ISM (Ferland 1996, and references therein), scaled by a factor of 1/2.5 to correct for sub-solar 
metallicity, were assumed. The resulting abundances relative to hydrogen are: 

He=1.00xl0"\ C=5.68xl0" 5 , N=3.73xl0" 5 , O=1.78xl0" 4 , Ne=4.68xl0" 5 , Mg=3.04xl0" 6 , Al=1.18xl0" 8 , 
Si=4.28xl0" 7 , S=6.48xl0" 6 , Ca=9.16xl0" 8 , and Fc=1.30xl0" 7 . 
Results for n# = 10 3 — 10 6 cm" 3 are given in models Al— A4 in Table 3. 

The observed [O III] A5007/H/3 flux ratio in the inner 20 pc core of the nucleus is matched with 
n#~10 4 — 10 5 cm" 3 . The [O I] A6300 luminosity can be generated if a high density component, n# ~10 6 cm" 3 , 
is also present (model A4). However, the other LINER features are not matched by these models due to the 
high ionization parameter. For example, the predicted [S II] A6726/H/3 ratio in model Bl is only 10% of the 
observed value and [N II] A6584/H/3 is underestimated by a factor at least two in all models. We find that 
no choice of n# or gas-phase abundances can reproduce these lines in the inner 20 pc region. 

Thus, if our derived luminosity and SED for the starburst is correct, much of the [S II] A6726 and 
[N II] A6584 emission must arise at further distances from the central starburst. In model B of Table 3, 
results are given for R=80, 130, and 180 pc. To satisfy the constraint imposed on the electron density by 
the [S II] A6716/[S II] A6731 ratio, we assumed n# = 10 3 cm" 3 . These results indicate R>130 pc is required 
to generate the [S II] A6726 and [N II] A6584 flux without overestimating the [O III] A5007 emission. Only 
at these large distances is the ionizing photon flux sufficiently dilute. However, this gas would lie near or 
beyond the edge of the area sampled by the aperture. Thus, any model of this extended emission region 
must account for aperture coverage effects. 

We have also computed photoionization models using the 6 Myr SED. Results are given in models C and 
D, which were computed with the same parameters as models A and B, respectively. Due to the decrease in 
photoionizing flux between the He II and He I edges compared to the 5 Myr SED (see Figure 5) , the 6 Myr 
model SED generates much weaker [O III] emission. The [O I] A6300 emission in the 6 Myr high density 
model (C4) is also much weaker because the SED is too steep to generate a sufficiently deep partially ionized 
zone in the backside of the cloud. For the same reason, the [S II] and [N II] emission in the extended emission 
line region (model D) is weaker than in the 5 Myr model. 

Best-Fit Model for the Emission Line Spectrum 

In Table 4, we present a composite model that represents our best-fit to the emission line spectrum of 
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NGC 4569 using the 5 Myr SED. A minimum of three distinct components is required for an acceptable fit. 
For component 1, we set n#=10 4,5 cm~ 3 at R=20 pc to reproduce the [O III] A5007/H/3 ratio measured in the 
core from the WFPC2 imagery. Component 2 consists of a low density gas at a relatively large distance from 
the ionizing starburst (low U) , which is required to produce the [S II] emission and much of the [N II] flux. 
We find that n#<900 cm~ 3 is needed to reproduce the [S II] A6716/[S II] A6731 flux ratio in the composite 
model. A high density component (component 3) is required to generate the [O I] A6300 luminosity. We 
chose component 3 to be at the location of the [S II] emitting gas, however, this component could be located 
closer to the starburst, as illustrated in model A4 in Table 3. The density for this component, n# = 10 5 cm~ 3 , 
was chosen to maximize [O I] A6300 and [N II] A6584 to provide the best overall fit. This composite model 
also satisfies the upper limits measured for the UV emission lines given in Table 1. Specifically, the predicted 
He II A1640 flux, listed in Table 4, is equals the upper limit. 

In deriving this model, the relative contribution of each component to the total emission was chosen to 
best match the observed spectrum. The covering factors for the three gas components were then derived by 
setting the sum of their H/3 luminosities equal to the observed value, Lnp =8xl0 39 ergs s _1 . We corrected for 
the aperture size in determining the total covering factors of the two extended components. Assuming that 
the emission from these components comes from a shell of emission at R=180 pc from the starburst, we find 
that the 160 pc x 320 pc aperture samples a total of 50% of this emission. Thus, the covering factors listed 
in Table 4 equal the H/3 luminosity derived for each component, divided by the H/3 luminosity predicted 
by the models, divided by the fraction of the total emission sampled by the aperture. This is only a rough 
correction for aperture coverage, since it assumes that a uniform gas surrounds the starburst and that the 
aperture is centered on this emitting shell. 

The high covering factor derived for our model indicates we may be observing the central starburst 
through the emission line gas. This may be responsible for the large gas column density in our line-of-sight 
to the starburst revealed by the strong interstellar absorption features and heavy extinction evident in the 
UV spectrum. 

Although the above model is a simplification, the three components represent physically distinct zones 
of emission. For example, the [S II] and [O III] components must be isolated and spatially distinct, with 
little line emitting gas between the two regions. Otherwise, the [O III] A5007/H/3 and [S II] A6726/H/3 ratios 
cannot be simultaneously fit. Furthermore, the large difference in densities required for the [O I] and [S II] 
regions indicates these lines arise in distinct clouds. 

5. Discussion 

5.1. Implications for the Starburst SED 

The results of our photoioinization modeling impose limits on the shape of the integrated SED of the 
central starburst in NGC 4569. From models A and B in Table 3, we find that the actual SED cannot be 
much harder than the 5 Myr model. Due to the high luminosity of the starburst and the low density of the 

5 II emitting gas, a harder SED cannot reproduce the S II/O III flux ratio and simultaneously satisfy the 
constraint on the location of the O III gas derived from the WFPC2 imagery - the ionization parameter is 
too high given the small aperture used to obtain the spectrum. 

This implies that models with very large populations of W-R stars, such as the solar metallicity 5 and 

6 Myr models generated with STARBURST99, are not correct for NGC 4569 since they have very hard 
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SEDs. This is consistent with our limit derived on the W-R population and metallicity in § 3.2. Similarly, 
stellar atmosphere models that predict stronger flux beyond the He I ionization edge for O and B stars than 
the Lejeune et al. (1997) models will produce a worse fit to NGC 4569. 

Conversely, the 6 Myr SED, with only a weak contribution from hot W-R stars, is too steep to generate 
the LINER features (see models C and D in Table 3). It lacks the extended tail of hard photons that 
is required to generate a warm partially ionized zone in the backside of the clouds where these lines are 
emitted strongly. The 6 Myr model represents a typical H I region SED. This indicates that W-R stars must 
contribute to some extent to the ionization in NGC 4569, as claimed by Barth & Shields (2000). 

5.2. Implications of Uncertainties on the Derived Model 

The uncertainty in the distance, D, to NGC 4569 introduces an uncertainty in the absolute values of the 
measured and derived parameters (e.g., luminosities and spatial extent of the emission line gas, R). However, 
since R oc D and L oc D 2 , the ionization parameter is independent of D (U oc L/R 2 ). For a given ionization 
parameter, photoionization calculations show that the emission line luminosities scale proportionally with 
the ionizing luminosity, except for models with very high U, which are not applicable to NGC 4569. Thus, 
a change in D does not affect the essential characteristics of our photoionization modeling results, namely, 
the derived densities, covering factors, and geometry of the emission line components. 

The largest uncertainty in our modeling results is due to the uncertainty in the extinction law, and 
thus ionizing luminosity, for the starburst in NGC 4569. We have tested various extinction curves from the 
literature. There are several observational constraints that place limits on the steepness of the extinction 
curve in NGC 4569. Specifically, the unreddened starburst cannot exceed the optical continuum luminosity 
we derived for the inner core region using WFPC2 imagery (§ 2.3) and the upper limit on the bolometric 
luminosity implied by the far-infrared data (§ 2.4). These luminosity constraints provide a lower bound on the 
steepness of the extinction curve. The Savage & Mathis (1979) curve, which we employed for our analysis, 
satisfies both of these constraints. However, we find that the Calzetti et al. (1994) extinction law, which was 
empirically derived for starburst galaxies, greatly overestimates these luminosities and is thus too "grey" to 
be applied to NGC 4569. Conversely, curves that are much steeper than the Savage & Mathis (1979) curve 
underestimate the ionizing luminosity required to generate the emission line luminosity, as evidenced by the 
high covering factor derived from our photoionization models. For example, the SMC curve of Hutchings 
(1982) results in a factor of two deficit in ionizing photons, and the very steep curve derived for the Seyfert 1 
galaxy NGC 3227 by Crenshaw et al. (2001) gives an even worse fit. Additionally, these steep curves would 
require that the emission line extinction, as determined from the Balmer decrement, is much larger than the 
extinction of the starburst. 



5.3. Physical Interpretation of the Model Results 

As demonstrated in § 4.2, a specific geometry for the nuclear gas is required if we are to accept stellar 
photoionization as the sole ionization and excitation mechanism in NGC 4569. A complete model, therefore, 
must provide a physical interpretation for the distinct spatial components. 

The mechanical energy due to frequent supernovae and strong stellar winds associated with the young, 
massive stellar population may create this structure. In this scenario, the [S II] emitting gas in NGC 4569, 
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component 2 in our best-fit model, may be swept-up interstellar gas surrounding the evacuated hot cavity 
of an interstellar supershell. The more compact [O III] A5007 clouds revealed in the WFPC2 narrowband 
imagery may be remnants of the original molecular cloud from which the starburst was born. In agreement 
with this, star-forming molecular clouds are thought to consist of dense condensations (e.g., n# ~10 2 - 
10 5 cm -3 ) embedded in very tenuous intervening gas (Efstathiou et al. 2000, and references therein). The 
diffuse component would be swept up relatively rapidly by the mechanical energy input by the young stars 
and may lead to the bi-component spatial structure proposed here. 

To assess the validity of this model, we calculated the initial ambient density that would be required to 
produce the line luminosity of the putative extended [S II] emitting region, assuming this gas is comprised 
of swept-up, diffuse material. Using L#/3=4xl0 39 erg s , which was derived for component 2 in our best-fit 
model, a minimum of 2xl0 61 hydrogen atoms are necessary to produce the emission. This is a lower limit to 
the total gas since it accounts only for the line emitting gas. Dividing this total by the volume enclosed in 
an R=180 pc region (consistent with our best-fit model), an average hydrogen density of n#=0.02 cm~ 3 is 
found for the tenuous component of the pre-starburst interstellar gas. This is typical of a diffuse interstellar 
density. 

The constraints derived on the location of the [S II] gas from our modeling results are consistent with 
the extent of the nuclear radio emission (see discussion in § 2.4). Thus, the extended radio flux could be 
free-free emission from the [S II] emitting nebular gas or non-thermal emission generated in shocks at the 
interface of the swept-up interstellar medium. Additionally, if we assume a typical temperature and density 
for the hot, tenuous cavity implied by this model, e.g., T=10 6 K and ny <10~ 2 cm~ 3 , we find that the 
X-ray luminosity predicted by the hot plasma models of Raymond & Smith (1976) is much weaker than the 
ROSAT measurement of the nucleus of NGC 4569. Thus, the X-ray emission is likely dominated by X-ray 
binaries and/or supernova remnants. 

Therefore, all observations are consistent with this model. We conclude that it provides a viable explana- 
tion for the bi-component gas required by our photoionization modeling, and that the stellar photoionization 
hypothesis remains a plausible solution for NGC 4569. 

6. Summary 

Our spectral population synthesis modeling results indicate the nuclear starburst in NGC 4569 is 5-6 Myr 
old and implies a relatively brief duration of star formation. The lack of detectable He II A1640 emission 
constrains the population of W-R stars and implies a sub-solar metallicity, with Z <Z Q /2.5. Luminosity 
constraints limit the steepness of the extinction curve that applies to NGC 4569. The Galactic curve of 
Savage & Mathis (1979) satisfies all luminosity constraints, and gives a similar reddening for the starburst 
and nebular line emission, E(B — V) =0.50. These results imply extreme conditions in the nuclear starburst, 
with ~5xl0 4 O and B stars compacted in the inner 9" x 13"region of the nucleus. Additionally, we find 
that this young, very dense starburst is distinct from the older, more spatially extended population of A 
supergiants detected in the optical spectrum (Keel 1996). 

The results of our photoionization modeling build on the previous analysis of Barth & Shields (2000), 
which demonstrated the general conditions required for a starburst to produce the enhanced LINER features 
observed in NGC 4569. Using all available observational constraints on the density, luminosity, and geometery 
of the nuclear region in our analysis, we have identified the physical elements that are required if the transition 
object spectrum in NGC 4569 is generated solely by stellar photoionization: 
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(a) At least two spatially distinct line emitting components are needed - a very compact region emitting 
[O III] A5007 and a more extended component for [S II] AA6716,6731. We propose a model in which these 
components correspond to the remnants of the parent molecular cloud of the starburst and to a swept-up 
"shell" of interstellar material, respectively. 

(b) A high density (n# ~T0 6 cm~ 3 ) gas is required to explain the [O I] A6300 emission. 

(c) The 5 Myr, Z=Z©/2.5 model SED derived with STARBURST99 (Leitherer et al. 1999) must be similar 
to the actual flux distribution of the starburst. If our derived extinction is correct, luminosity constraints 
limit the SED from being much harder. Conversely an SED that is much softer cannot generate the 
LINER features. These results are consistent with our population synthesis results, which indicate Zq/IOk 
Z <Z /2.5. Thus, our photoionization models support the assertion by Barth & Shields (2000) that W-R 
stars play an important role in the ionization of NGC 4569. 

In conclusion, we have shown through detailed modeling that photoionization by the nuclear starburst 
is capable of generating the observed LINER/H II region transition object spectrum, consistent with all 
observations. Our model places specific requirements on the physical conditions in the nucleus. Thus, it 
makes several predictions that can be tested with future observations. 

Facilities of the Laboratory for Astronomy and Solar Physics at NASA/GSFC were used in performing 
this research. We thank the referee, W. Schmutz, for providing very helpful comments for this paper. We 
also thank Steve Kraemer and Mike Crenshaw for useful discussions and advice, Margaret Smith Neubig 
for assisting with the population synthesis analysis, and Gary Ferland (CLOUDY) and Claus Leitherer 
(STARBURST99) for making their software publicly available. 
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Fig. 1— HST/WFPC2 continuum image of NGC 4569 using the F547M filter. The square delineates a 
0".35 x 0".35 region corresponding approximately to the spatial extent of the compact starburst. For the 
adopted distance to NGC 4569, each side of the square represents about 25 pc. 

Fig. 2. — Multi-wavelength spectrum of the nucleus of NGC 4569. The 6 cm and 2 cm VLA radio data (*) 
are upper limits placed on any unresolved (1") core emission from Neff & Hutchings (1992). The IRAS 
data points (o) are from Neff & Hutchings and sample the inner The UV (dot-dashed line) is from 

the HST/FOS spectrum. Barth et al. (1998) found the UV source has a FWHM size of w0".15. Optical 
data (solid line) are for the inner 0".35 core of emission measured from the F547M WFPC2 image (this 
study). The X-Ray data (dashed line) is a power law fit to the observed ROSAT/RRl data, which has «4" 
resolution, and thus represents an upper limit to the inner core luminosity. Upper limits are indicated with 
arrows. 

Fig. 3. — Population synthesis spectra for various ages plotted for the (a) N V A1240, (b) Si IV A1400, (c) 
C IV A1550 and (d) He II A1640 spectral regions. Both the theoretical (solid line) and observed spectrum 
(dashed line) have been normalized. Empirical fits to the broad wind absorption features are plotted on the 
spectra for Si IV and C IV (smooth line). All model spectra were generated with CLUST_FLUX using £ms 
and Zq. The prominent stellar wind and photospheric features are identified above the spectrum. Strong 
interstellar lines are marked (X) below the spectrum. The spectrum of NGC 4569 has been convolved 
(smoothed) with a a~ 1 A Gaussian filter. All model spectra were smoothed with a er=1.5 A filter to best 
match the resolution of the NGC 4569 data presented in the figure. 

Fig. 4. — Population synthesis spectra for 1890-1990 AThe model spectra (solid line) and the observed 
spectrum (dashed line) have been normalized. The Fe III lines are identified at the top of the and 5 Myr 
models with tick marks. The spectra were smoothed as described in the caption of Figure 3. 

Fig. 5. — 5 Myr (solid line) and 6 Myr (dashed line) model starburst spectral energy distributions. The 
SEDs were derived with Z=Z Q /2.5 and the Miller & Scalo (1979) IMF, assuming an instantaneous burst of 
star-formation using the population synthesis code STARBURST99 (Leitherer et al. 1999). 
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Table 1. 


Emission Line 


Data 


Line 


Observed 3 - 


Dereddencd b 


H/3A4861 


1.0 C 


1.0 d 


[0 III] A5007 


1.2 (± 0.2) 


1.1 (± 0.5) 


[0 I] A6300 


0.31 (± 0.10) 


0.19 (± 0.08) 


Ha A6563 


5.0 (± 1.3) 


2.9 (± 1.1) 


[N II] A6583 


4.5 (± 1.0) 


2.6 (± 0.9) 


[S II] AA6716,6726 


2.0 (± 0.4) 


1.2 (± 0.4) 


C IV AA1548,1551 


<0.04 


<0.3 


He II A1640 


<0.04 


<0.3 


C III] A1909 


<0.02 


<0.2 


C II] A2326 


<0.02 


<0.2 


[0 II] A2470 


<0.02 


<0.1 



a Optical fluxes from Ho et al. (1997a); UV upper 
limits from Maoz et al. (1998); Fluxes relative to H/3 

b Corrected for E(B-V)=0.48, Savage & Alathis 
(1979) extinction curve 

c Observed: L^^^l.exlO 39 erg s _1 ; assuming D = 
16.8 Mpc 

d Dereddcned: L^f^^S.lxlO 39 erg s _1 ; assuming D 
= 16.8 Alpc 
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Table 2. FOS Observations of NGC 4569 



Data Set Grating Aperture Resolution Exposure 

A sec 



Y3400104T G130H 0.86 -1.0 2370 

Y3400105T G130H 0.86 -1.0 2420 

Y3400106T G190H 0.86 -1.5 1320 

Y3400107T G270H 0.86 -2.0 690 



- 23 - 



Tabic 3. Photoionization Model Results 



Model 


SB Ago a 
Myr 


R 

pc 


nff 

cm~ 3 


Log(U) 


10 40 org s- 1 


[O III] A5007 b 


[O I] A6300 


[N II] A6583 


[S II] A6726 c 


R[s//] d 


Al 


5 


20 


10 3 


-1.5 


0.97 


6.6 




0.19 


0.11 


0.84 


A2 






10 4 


-2.5 


1.7 


4.6 


0.01 


0.95 


0.25 


0.70 


A3 






10 5 


-3.5 


2.0 


0.78 


0.12 


1.2 


0.24 


0.50 


A4 






10 6 


-4.5 


2.1 




0.47 


0.38 


0.10 


0.40 


Bl 


5 


80 


10 3 


-2.7 


1.0 


3.5 


0.02 


1.3 


0.77 


0.84 


B2 




130 


10 3 


-3.1 


2.1 


2.0 


0.04 


1.8 


1.3 


0.84 


B3 




180 


10 3 


-3.4 


2.1 


0.83 


0.09 


2.1 


1.8 


0.84 


CI 


6 


20 


10 3 


-1.4 


1.8 


2.0 




0.79 


0.10 


0.84 


C2 






10 4 


-2.4 


2.1 


0.77 


0.01 


1.4 


0.24 


0.70 


C3 






10 s 


-3.4 


2.5 


0.07 


0.04 


0.93 


0.13 


0.50 


C4 






10 6 


-4.4 


2.6 




0.16 


0.18 


0.03 


0.40 


Dl 


6 


80 


10 3 


-2.6 


2.3 


0.43 


0.01 


1.4 


0.65 


0.84 


D2 




130 


10 3 


-3.1 


2.4 


0.15 


0.02 


1.5 


1.0 


0.84 


D3 




180 


10 3 


-3.3 


2.4 


0.07 


0.03 


1.5 


1.3 


0.84 



"Model starburst SED from Leither et al. (1999) 
b All fluxes relative to H/3 
"Total [S II] A6716+A6731 flux 
R[S//]=L[ S7 j] A6716 /L[ S/J -] A6731 
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Tablc 4. Best-Fit Model using the 5 Myr Starburst SED 





R 
pc 


riff 
cm~ 3 


Log(U) 


[O III] A5007" 


[O I] A6300 


[N II] A6583 


[S II] A6726 b 


He II A1640 


R[sff] c 


Covering Factor 


Comp 1 


20 


1Q 4.5 


-3.0 


2.9 


0.04 


1.3 


0.26 


0.41 


0.5 


0.14 


Comp 2 


180 


900 


-3.4 


0.83 


0.09 


2.1 


1.8 


0.30 


0.84 


0.38 


Comp 3 


180 


10 5 


-5.4 




1.1 


1.1 


1.6 


0.01 


0.54 


0.13 


Composite 








1.3 


0.24 


1.8 


1.3 


0.30 


0.80 




Observed 








1.1(±0.5) 


0.19(±0.08) 


2.6(±0.9) 


1.2(±0.4) 


<0.3 


1.0(±0.2) 





a All fluxes relative to H/3 
b Total [S II] A6716+A6731 flux 

°R[Sff] =L[Sf f ] A6716/L[Sf f]A6731 

d Lff comprised of 33% Comp 1, 50% Comp 2, 17% Comp 3 
e Dereddened by E(B-V)=0.48 using Savage & Mathis curve 



